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1. INTRODUCTION 

Recently, the computer aided diagnosis (CAD) is widely spread which it enhances the understanding 
of the medical images. Edge detection (ED) plays the important roles in the field of image processing 
particularly in the segmentation and clustering procedures. It is based on calculating the changes of the 
intensity for the neighbour pixels in the image. ED has been widely used in the field of medical image 
processing [1]-[4]. In the paper, Anas et al. studied the performance of several edge detection methods on 
different X-ray images [1]. They stated that the best Gaussian noise removing algorithm in X-ray medical 
images was that median filter. Study by Chucherd introduced a generalized gradient vector flow method as 
an edge detection procedure to detect a breast cancer in ultrasound images [2]. Chucherd mentioned that the 
proposed method produced higher segmentation accuracy comparing with as Sobel and Canny segmentation 
methods. The neural network with the aid of conventional edge detection method was also used for 
segmentation blood vessels of retina [3]. Fuzzy geometrical approach has also been used to enhance the 
contrast of digital images for better classification [5]. In this study, the conventional edge detection filters 


Journal homepage: http://ijeecs.iaescore.com 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 m) 87 


such as Sobel and canny filters were used to extract the neural network algorithm features. Most of the state 
of art studies used the Sobel operator as a tool for edge detection and image segmentation [2], [3]. 

The Sobel filter is one of the most effective and popular methods in ED and image segmentation [6]. 
It is based to detect the vertical and horizontal edges in the image. There are several platforms, programs and 
languages that can be used to implement the Sobel operator such as MATLAB Software [7], field 
programmable gate array (FPGA) [7], [8], and OpenCV based FPGA [9], [10]. In terms of FPGA, there are 
various FPGA studies have implemented the Sobel edge-detection using different algorithms. Nausheen et al. 
proposed an implementation of the Sobel filter using very high speed integrated circuit (VHSIC) hardware 
description language (VHDL) [11], [12]. Moreover, Zhang et al. designed a Sobel operator using VHDL 
[13]. Menaka et al. also implemented an ED system based on FPGA [14]. All these mentioned studies aimed 
to faster the performing processes and to reduce the spatial complexity of the FPGA [11], [13], [14]. Sobel 
filter has been also used by Solanki and Tailor to design a video streaming architecture using high level 
synthesis [15]. In Solanki and Tailor study, zybo board (based also on ZYNQ 7000) was used, which has 
provided adequate peripherals for implementation. Although, the Sobel filter has been used in most of the 
edge detection studies, these studies used the standard codes for implementation of the Sobel filter. This 
study aims to parallel implementation of the Sobel operator based on the customization code of FPGA and to 
compare the results in terms of quality of edge detection, spatial complexity and time consuming with state- 
of-art studies. 

The organization of this research is as follows: methods and material which highlights the 
theoretical aspect of Sobel filter mathematics. It also depicts the platform used in implementing of the 
hardware part of the system. Practical system section includes the practical part for implementing the 
proposed techniques with their hardware architectures that are used to detect the edges of the cells in blood 
images. Result and discussion section provides empirical results and performance analysis obtained by 
implementing methods used in the proceeding sections. Finally, conclusion section presents the conclusions 
of the study. 


2. METHODS AND MATERIAL 

The first part in this section mainly presents the edge detection techniques used in this paper. 
Generally, edge detection techniques are based on the concept of 2D spatial filtering operation. Spatial 
filtering operation refers to periodicity with which pixel values change. Just contrary to other techniques of 
the edge detection, the algorithms of first difference-based operator are the simplest to implement on 
hardware. Based on that, starting from the upper left side of the image, horizontal and vertical masks slide all 
over the image. Also, the masks move in the right direction on the image until they arrive at the end of each 
row. Then, they start again from the left most next row of the image [16]. To include all image pixels by the 
filter kernel, it is to a usual procedure to pad the image by the number of rows and columns depending on 
kernel dimensions. This had not been followed here because the main information dose not reside in the 
image edges but anywhere else. The second part in this section shows the materials that are used in 
implementation of this filter to detect the edges of the blood cells images. 


2.1. Sobel edge detector 

Sobel operator is a classical first order edge detection operator. It is used to compute the 
approximation of the image intensity function. At every point in the image, the result of the Sobel operator is 
the corresponding norm of this gradient vector. It makes use of two 3x3 masks convolved with the original 
image and calculates the image approximation. Both convolution masks are designed in such a way to detect 
the edge in both horizontal and vertical directions. The following (1) and (2) are used to calculate the G, and 
G, masks (the gradient component of each orientation G, and Gy) [16]. 
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The separate output gradients is combining in order to obtain the absolute magnitude of the gradient 
at every points of the image. In (3) illustrates the gradient magnitude |G|. The direction of the gradient is 
calculated by (4), while Figure 1 illustrates the flowchart of the steps followed in the design of this filters. 
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IG| = [G.2 + G,? (3) 
6 =arctan (5) (4) 


Finally, binary images can be created by changing the threshold value. The threshold value is 
obtained by replacing each pixel with either a black or white one. The replacing of black pixel is dnoe when 
the image intensity does not exceed some fixed constant T (that is, I (i, j) < T), while white pixel whose 
intensity exceeds that constant. 
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Figure 1. Flowchart of Sobel filter algorithm 


2.2. Implementation platforms 
The following platforms are used in this study: 

— MATLAB: the first software platform used in this work is matrix laboratory (MATLAB). The image 
processing toolbox (IPT) is being the most relevant to this study [17]. MATLAB (version 2019) is used 
to design Sobel filter which is used to manipulate image processing and to produce edge detection for the 
original image. 

— Xilinx OpenCV: this software platform is a template library optimized for FPGA high-level synthesis 
(HLS) which is used to create the image processing pipelines simply because it contains common 
algorithms such as: color conversion, image resizing, border and edge detection algorithms (Canny, 
Sobel), warp transformation, hog transform, matrix-matrix operations and others. It is completely 
possible to integrate Xilinx OpenCV to the Vivado HLS tool for accelerator-oriented projects. This 
library is used to build the first architecture that implement the Sobel edge detector [18]-{20]. 

— Vivado Design Suite: it is a software suite produced by Xilinx for synthesis and analysis 
of HDL designs. Xilinx Vivado 2018.3 is used to design the second architecture for detecting the edge of 
the image using Sobel filter [20]. 

— FPGA: The ZC702 board or Zynq-7000 XC7Z020 AP SoC is used to implement the proposed design of 
this study as a hardware platform. It comprises of a processing system (PS) organized around a dual-core 
ARM Cortex-A9 processor, and programmable logic (PL), that equivalent of traditional FPGA with 
additional new features such as integrated memory, variety of peripherals, and high-speed 
communications interfaces [21]-[24]. Both architectures designed (the first one using Xilinx OpenCV 
and the second one using Vivado Design Suite) are implemented using this hardware platform. 


2.3. Proposed system 


This section is divided into three parts: the first part deals with design the Sobel filter to detect edges 
using the MATLAB program, this technique is used for the purpose of comparing the results with the second 
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and third methods. The second part presents our first architecture design using the OpenCV library from 
Xilinx, where the build-in functions are used in the design of Sobel filter. The third part, where the Sobel 
filter algorithm has been programmed using HLS to produce a new special architecture built to detect edges 
in blood cells of microscopic blood samples images. 


2.3.1. Matlab based Sobel edge detection system 

In order to evaluate the proposed hardware system, function simulations have been done using 
MATLAB package and the results have been examined before the designs have been implemented on the real 
device. The Sobel filter algorithm is programmed manually step by step without using build-in functions, this 
is done for fair comparison with the proposed architectures in next sections. 


2.3.2. Xilinx OpenCV based Sobel edge detection system 

An architecture is built to detect the edges of cells in blood microscopic images. The architecture is 
creaed by using OpenCV library built-in functions. Figure 2 illustrates the flowchart of the Sobel filter 
implementation process using OpenCV. 
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Figure 2. Flowchart of the Sobel filter implementation process using OpenCV 


At this point, it is worth noting that the final gradient calculation after each step of the convolution 
process (G in (3) using square root) is done using the special addition function of this library (named 
AddWeighted) where different weights are given to images so that it gives a feeling of blending. Figure 3 
shows the block diagram obtained from implementing of this algorithm on the ZYNQ hardware platform. As 
shown in this Figure 3 the edge detection system designed using OpenCV consists of the Sobel IP core which 
is designed inside PL to detect the edges of an input source image and the ZYNQ processor inside PS which 
is used to connect this core through automated X-ray inspection (AXI) bus to input output peripherals. 
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Figure 3. Implementation of the Sobel algorithm on the Zynq chip using OpenCV 


2.3.3. Customized parallel HLS based Sobel edge detection system 

This part of the research presents the customized parallel HLS algorithm. Figure 4 shows the 
Flowchart of the Sobel filter parallel implementation process using HLS. It can be noticed that the code is 
written to work in parallel to find the gradient for each direction of the image. The gradient component of 
each orientation G, and Gy are calculated simultaneously then final G is calculated to reach our goal for 
finding the edges of the image faster. The parallelism process is illustrated in Figure 5, where the designed 
architecture consists of two IP core (named Sobel-Ix and Sobel-Iy) of the side PL. The other blocks of this 
design complete the work as in previous design; ZYNQ processor of PS, AXI buses and control. 
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Figure 4. Flowchart of the Sobel filter parallel implementation process using HLS 
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Figure 5. Implementation of the Sobel algorithm on the Zynq chip using Customized parallel HLS 


3. RESULTS AND DISCUSSION 
3.1. MATLAB based Sobel edge detection system results 

Table 1 represents the results obtained using MATLAB for microscopic blood images. Where the 
upper images represent the original input image. The lower images used in this table are represent the 
resulted edge detected images using MATLAB method. 


Table 1. MATLAB based Sobel Edge detection system results 
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3.2. Xilinx OpenCV based Sobel edge detection system results 

Table 2 represents the results that were obtained using OpenCV library for microscopic blood 
images. Where the upper images represent the original input image. Images in the lower part of the table 
represent the resulted edge detected images using Xilinx OpenCV method. 

Figure 6 shows the overall system hardware utilization resources that used to implement the hardware 
Sobel Edge detection system. The examining resources are based on Zynq-7000 XC7Z020 board. It can be 
noticed that the maximum utilization is 18% go for the DSP block which is used to for convolution process. 


Table 2. Xilinx OpenCV based Sobel Edge detection system results 
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Resource Utilization Available Utilization % 
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LUTRAM 172 17400 0.99 
FF 7411 106400 6.97 
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Figure 6. OpenCV based Sobel Edge detection system hardware utilization resources 


3.3. Customized parallel HLS based Sobel edge detection system results 

Table 3 represents the results obtained using customized parallel HLS algorithm. Where the upper 
images represent the original input image, while the lower images show the resulted edge detected images 
using customized parallel HLS method. Figure 7 shows the customized parallel HLS based Sobel Edge 
detection system hardware utilization resources of Zynq-7000 XC7Z020 board that are used to implement 
this design. By examining this figure, we can see that the maximum utilization is 14% go for the LUTs. 


Table 3. The ccustomized parallel HLS based Sobel Edge detection system results 
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Figure 7. Customized parallel HLS based Sobel Edge detection system hardware utilization resources 
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3.4. Performance evaluation results 

The performance evaluation of the proposed system has been performed by calculating the mean 
square error (MSE) and peak signal-to-noise ratio [25], [26]. According to the obtained mean square error 
(MSE) results demonstrated in Table 4, all techniques demonstrate very similar results, but MSE obtained 
from the OpenCV based method was better performance than the others. The same conclusion can be noticed 
in the Table 5 which represents the peak signal-to-noise ratio in decibels (dB). Finally, Figures 8 and 9 show 
the estimated time of the core using OpenCV and customized parallel HLS, respectively. It can be noticed 
that the built-in functions algorithm has higher estimated time, which is critical and leads to a delay of 
obtaining the results, compare to the estimated time of the proposed customized algorithm with no or smaller 
delay. 


Table 4. The mean square error (MSE)x 10+ 
Image No. 1 2 3 4 5 6 7 8 
Technique 

Matlab 37 29 39 3.7 36 34 37 3.7 
OpenCV 36 27 38 36 36 34 3.7 3.7 
Customized 3.9 29 40 38 38 36 39 3.9 


Table 5. The peak signal-to-noise ratio in decibels (dB) 


Image No. 1 2 3 4 5 6 7 8 
Technique 
Matlab 24 34 21 24 25 2.7 23 24 


OpenCV 20 3.2 2.2 25 25 2.7 24 25 
Customized 2.2 34 20 22 23 24 21 2.1 


=| Timing (ns) 


Clock | Target | Estimated | Uncertainty 


ap_cik | 10.00 | 10.049 | 1.25 


= Latency (clock cycles) 


= Summary 
| Latency Interval | 
min max i min max | Type 
| 134699 | 141225 | 134688 | 141209 | dataflow 


Figure 8. OpenCV based Sobel Edge detection system performance estimation 
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Figure 9. Customized parallel HLS based Sobel Edge detection system performance estimation 


4. CONCLUSION 

The microscopic blood images play the important roles to diagnose different biological diseases. 
However, the diagnostic error can be produced due to human error. Therefore, this study aimed to implement 
an automated edge detection system for microscopic blood images. It proposed a new method of parallel 
implementation for Sobel operator based on the customization code of FPGA. It is observed that the FPGA 
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customized codes give the best visual and less timing results than the ready to used algorithms such as 
MATLAB, using OpenCV standard code. The proposed system can not only apply to the microscopic blood 
images but also to other medical images, such as X-ray, computer tomography and magnetic resonance 
images. 
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